Extracting quantum work statistics and fluctuation theorems by single qubit interferometry 
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We propose an experimental scheme to verify the quantum non- equilibrium fluctuation relations using current 
technology. Specifically, we show that the characteristic function of the work distribution for a non-equilibrium 
quench of a general quantum system can be extracted from Ramsey interferometry of a single probe qubit. 
Our scheme paves the way for the full characterisation of non-equilibrium processes in a variety of quantum 
systems ranging from single particles to many-body atomic systems and spin chains. We demonstrate our idea 
using a time-dependent quench of the motional state of a trapped ion, where the internal pseudo-spin provides a 
convenient probe qubit. 

PACS numbers: 03.67.Mn, 03.67.Lx 



Introduction - In microscopic systems, non-equilibrium be- 
haviour is dominated by thermal and quantum fluctuations. 
Obtaining a comprehensive understanding of these phenom- 
ena is therefore of great fundamental importance. In the past 
two decades, the discovery of the non-equilibrium fluctuation 
relations has made a significant contribution to this endeavour 
by characterising the full non-linear response of a microscopic 
system subject to a time-dependent force J3-0] . Initially, these 
relations were derived for classical systems and experimen- 
tally confirmed in single-molecule stretching experiments 0]. 
More recently, their extension to the quantum regime has lead 
to theoretical progress in understanding the microscopic un- 
derpinnings of the laws of thermodynamics |@, 0]. Notably, 
however, an experimental verification of the quantum fluctua- 
tion relations is still forthcoming. 

In this letter, we propose an experimental scheme to extract 
the full statistics of work done in a non-equilibrium transfor- 
mation of an arbitrary closed quantum system. The crux of 
our proposal is that the characteristic function of the work 
distribution can be extracted via Ramsey interferometry of a 
suitable probe qubit. Extracting the statistics of work in this 
way circumvents the requirement to perform projective energy 
measurements on the system of interest, as in previous pro- 
posals (Hi, making our scheme generally applicable to a range 
of systems, including Bose if^L Tioll and Fermi gases 111, 121. 
spin chains [13] and quenched ion strings [14]. A further 
advantage of our scheme is that the temperature of the sys- 
tem is extracted directly from the measurement signal. This 
feature suggests potential future applications for in-situ, non- 
destructive thermometry in systems such as Bose-Einstein 
condensates. We demonstrate the feasibility of our proposal 
using realistic parameters for the example of a trapped ion 
interacting with an external laser field. In this case, the inter- 
nal pseudo-spin state of the ion provides a convenient probe 
qubit (see Ref. Il5ll for details of a similar scheme using hy- 
brid opto-/electro-mechanical devices). 

Non- equilibrium quantum thermodynamics - We begin by 
briefly reviewing some important concepts from non- 



equilibrium thermodynamics and defining the quantities that 
will be of interest in the rest of this letter. We therefore 
consider a closed quantum system described by a Hamil- 
tonian H(X) containing an externally controlled parameter 
A(t). At time t = the control parameter has the initial 
value A(0) = A; and the system is prepared in the Gibbs 
state Qp(\) = exp[—0H(\i)]/Zp(\i), where Zp(X) ■= 
tr[exp(— /3H (A))] is the partition function at inverse temper- 
ature /3. The system is driven away from equilibrium by vary- 
ing A(t) in a pre-defined, but otherwise arbitrary way over the 
quench time interval £q to its final value A(4q) = Af. The ini- 
tial and final Hamiltonians have the spectral decompositions 
H(k) = En e n\n){n\ and H(X f ) = £ m e m \m)(m\, respec- 
tively, and the protocol H(Xi) — > H(Xf) that connects them 
generates the unitary evolution E/(£q). 

The work done on the system W is defined by two projec- 
tive energy measurements: The first, at t = 0, projects onto 
the eigenbasis of H(Xi) and gives the outcome e n with a prob- 
ability p n = exp(— /3e n )/2 | g(Ai). The second measurement, 
at the end of the quench protocol t = iQ, projects onto the 
eigenbasis of H(Xf) and gives the outcome l m with probabil- 
ity p m \n = I (^{U (tQ)\n) | 2 . The quantum work distribution, 
which encodes the random fluctuations in non-equilibrium 
work arising from both thermal (p n ) and quantum measure- 
ment (p m \n) statistics over many identical realisations of the 
quench protocol, is thus, 

P F (W) := ^2p n p m \n${W - (e m - e„)). 

Here 'F' denotes that this is the work distribution for the for- 
ward process H(X\) — > -ff (Af). The corresponding backward 
work distribution Pb (W) is obtained by preparing the system 
in the Gibbs state Qp(Xf) of the final Hamiltonian and subject- 
ing it to the time-reversed protocol H(Xf) — > H{X\) generated 
by the evolution U'(to). 

By studying the fluctuations in non-equilibrium work it is 
possible to extract important equilibrium information. This is 
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revealed by the non-equilibrium fluctuation relations, in par- 
ticular, the Tasaki-Crooks relation 0] 



Pf(W) 
Pb(-W) 



a P(W-AF) 



(1) 



which shows that, for any closed quantum system undergoing 
an arbitrary non-equilibrium transformation, the fluctuations 
in work are related to the equilibrium free energy difference 
AF = (l//3)ln[Z /3 (A i )/Z /3 (Af)]. This relationship is further 
emphasized by a corollary to Eq. ([1} known as the Jarzynski 
equality J2l; 
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which states that AF can also be extracted from the properties 
of the forward (or backward) work distribution alone. 

The primary quantities of interest in this letter are the 
forward and backward characteristic functions, defined as 
the Fourier transform of their corresponding work distribu- 
tion 11611 . Hence, the forward characteristic function (taking 
h=l) 

X f(u) := J dWe iuW P F (W), 

= tr[W(t Q )e iu6 ^U(t Q )e- iu6 ^g (X{)], (2) 

while the backward characteristic function Xb(u) '= 
J dWe mW P B (W). 

Previous experimental proposals to extract the full statistics 
of work, and hence verify the quantum fluctuation theorems, 
have sought to directly measure the work distribution via a 
series of projective energy measurements JH. However, even 
for quantum systems of modest complexity this can be practi- 
cally challenging. In the following section, we show how this 
difficulty can be avoided by instead extracting the character- 
istic function of the work distribution using well-established 
experimental techniques. 

Experimental extraction of the characteristic function - The 
purpose of our proposal is to measure the work done in a non- 
equilibrium transformation of a generic quantum system by 
temporarily coupling it to an easily-addressable probe qubit. 
We assume that the total Hamiltonian describing the probe 
qubit and system of interest has the form Hj(t) = —fr z + 
Hs + H\(t) where A is the splitting between the ground \l) 
and excited |^) states of the qubit, which are eigenstates of 
the spin- 1/2 Pauli-z operator a z (similarly a x and fr y denote 
the Pauli-x and -y operators) and H$ is the time-independent 
bare Hamiltonian of the system of interest. The qubit-system 
interaction term Hi(t) contains all of the time-dependence and 
is assumed to have the form 

Hi(t) = U(t)\l)(l\ + <n(t)tt)(t\) ®v, 



where g^(t) and g±(t) are externally controlled parameters 
and V is a perturbation acting on the system of interest. 



With an interaction of this form, an isolated quench of the 
system described by the protocol H(Xi) — i?s + \V — ► 
H(\f) = Hs + XfV can be realised by varying both of the 
spin-dependent parameters, g±(t) and <?-(-(£), from the initial 
values .94.(0) = g T (0) = A; at t = to the final values 
9l(t(l) = dti^o) = Af at £ = £q according to the same, pre- 
defined protocol. 

To extract the characteristic function for this quench, how- 
ever, a related but distinct extraction protocol must be per- 
formed on the system and qubit as part of a modified Ram- 
sey sequence in which gi(t) and gf(£) are varied over the 
Ramsey time interval £r. Crucially, during the extraction pro- 
tocol, the parameters 174. (£) and gf (£) will differ so that the 
probe qubit and system become coupled. Explicitly, the ex- 
perimental procedure to extract the characteristic function is 
as follows: i) For times t < the qubit is decoupled from 
the system by holding the spin-dependent couplings fixed at 
gi(0) = 3^(0) = A;. Furthermore, the qubit and system are 
thermalised in the product state p = \l) (\.\ <g) gp(\i) by ensur- 
ing that j3A S> 1. ii) At £ = 0, a Hadamard operation <7h = 
(fr x + a z )/\/2 is applied to the qubit. Hi) The spin-dependent 
couplings are independently varied over the Ramsey time in- 
ternal £r from their initial value 94,(0) = <?t(0) = A; to the 
final value g± (£ R ) = .94- (£r) = Af. This protocol generates the 
unitary evolution operator T(£r) that acts in the joint Hilbert 
space of the qubit and system to generate a conditional dy- 
namical quench of the system contingent upon the state of the 
probe qubit. For the sake of clarity, we stress that this pro- 
tocol is not the same protocol for which we wish the extract 
the work done, iv) At the end of the protocol, t = £r, the 
qubit and system are automatically decoupled and a second 
Hadamard operation is applied to the qubit at. 

The output state of the probe qubit at the end of the Ramsey 
sequence is, thus, 
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where we have introduced the decoherence factor 

L(£ R ) =tr s [T t t (£ R )T 4 (£ R )^(A i )]. 



ItXtl, (3) 



(4) 



Here, the unitary operators Tj,(£r) = (||T(£r)||) and 

T^{t R ) = (t|T(£ R )|t) act in the Hilbert space of the sys- 
tem of interest and describe its evolution under the two dif- 
ferent time-dependent quenches generated by <?4.(£) and g^(t), 
respectively. Consequently, the Ramsey sequence (shown in 
Fig. [TJ a )) creates an entangled state between the basis states 
of the probe qubit and the two quenched states of the sys- 
tem, f^gpiX^fl and f t [^(Ai)]f|. The real U[L(t K )] and 
imaginary 9[X(£ R )] parts of the decoherence factor define the 
populations and coherences of the probe qubit density matrix 
in Eq. (01 and are experimentally reconstructed by measuring 
a z and fry over many identical experimental runs. 
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A direct relationship between L(£r) in the Ramsey scheme 
and the characteristic function xf( u ) f° r the quench protocol 
H(\i) = H s + AjV -> i?(A f ) = H s + A f y is established 
by judiciously engineering the unitary operators T±(t R ) and 
2\-(4r); As illustrated in Fig.[TJb)-(c), the spin-dependent con- 
trol parameters g\(t) and g+(t) are varied so that the system 
state entangled with |t) undergoes the quench A; — > Af over 
the time interval £q, followed by constant evolution at Af up to 
the Ramsey time £r. Simultaneously, the system state entan- 
gled with |t) undergoes constant evolution at Aj followed by 
the quench. The corresponding unitary operators are thus 



(a) 



I) — <*h 



2t(tn) 



-i(t R -t Q )H(\,) 



U(t Q ), 



which, after identifying u = t R — tQ, show that the deco- 
herence factor in Eq. (@) coincides exactly with the forward 
characteristic function in Eq. (|2). Hence, the characteris- 
tic function is extracted by embedding the evolution U(to) 
into the qubit-system evolution and repeating the protocol for 
different run times £r > tQ. The corresponding backwards 
characteristic function is obtained by a straightforward mod- 
ification of the above scheme. In both cases, the work dis- 
tributions Pp(W) and Pb(W) are obtained from the inverse 
Fourier transform of their respective characteristic functions 
(cf. Eq. (f2]i). In the following section we illustrate how this 
scheme can be implemented using a conventional trapped ion 
arrangement under realistic circumstances. 

Implementation using a trapped ion - We consider a sin- 
gle ion of mass M contained in a linear Paul trap 111 811 . By 
using the S\/2 ground state Zeeman sublevels of the ion 
|m = 1/2) = |t) and |m = -1/2) = ||) this system pro- 
vides an ideal realization of a spin- 1/2 particle confined in a 
harmonic potential. We therefore have Hs = uJo(a'a + 1/2), 
where ujq is the natural frequency of the oscillator and a' (a) 
is the oscillator raising (lowering) operator. The trapped ion 
setup has a number of distinguishing features: First, accurate 
detection of the ion's internal states can be accomplished by 
observing the scattered fluorescence from near-resonant driv- 
ing of a cycling transition. Second, transformations between 
internal states can be implemented by a Raman transition (e.g. 
performing the Hadamard operation cth via a ir/2 pulse) and 
the tunable azimuthal phase of the transition permits both 
(<t z ) and (a y ) to be determined from the fixed final measure- 
ment 01811 . Third, precise preparation of the initial thermal 
state Qp, with mean phonon number n — [cxp(/3o;o) — 1] — 1 , 
can be achieved by allowing heating after resolved-sideband 
laser cooling to the motional ground state or Doppler cooling 
on the Si/2 to P\ji transition llj, 181. 

Similar to an earlier proposal studying quantum chaos us- 
ing a trapped ions II 1711 . we quench the motional state of the 
ion by illuminating it with a far-detuned elliptically polarized 
standing wave laser field. Further, since the a + and o~ polar- 
ized contributions couple exclusively to the |t) and |f ) states, 
respectively, they induce a spin-dependent optical dipole po- 
tential for the ion 111 911 . Hence, after making the rotating- wave 
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FIG. 1 : (a) The Ramsey sequence represented as a quantum circuit. 
The probe qubit in the upper branch is prepared in the \\) state and 
the system of interest is prepared in the state Qp(Ai) defined in the 
text. A black (white) circle indicates that the operation is controlled 
on the probe qubit being in the |t) (|t)) state. In (b) and (c) examples 
are shown of the time variation of the spin couplings gi(t) and Qf(t) 
over the Ramsey scheme time tn = tg + u required to obtain the 
characteristic function xf( u )- 0>) A forward process of the form 
X(t) = Xi + (A/ - Ai) [l + tanh(t/T)]/2, where T is the switching 
time and the total quench time is tQ — 8T is shown, (c) A forward 
process is shown which is composed of repeated fast and slow tanh 
switching between \i and \ f, like that in (b), before ending at A/. 



approximation and adiabatically eliminating the far-detuned 
excited states, we find the interaction Hamiltonian 112811 



9l(t)\i)(i\ + StWItXtl ® shAH + 0), 



where k is the magnitude of the wave-vector orientated along 
the axis of the trap for both polarisations, and <f> is the phase 
of the standing waves relative to the trap centre at x = 0. 
The couplings f2-|-(i) and fli(t) are the time-dependent Rabi 
frequencies, which are independently controlled by varying 
the laser intensity for the corresponding polarization. 

In the Lamb-Dicke regime, quantified by r\ = kxo <C 1 
where xq = (2Mujo)^ 1 ^ 2 , the extent of the ions motion is 
small compared to the spatial variation of the optical dipole 
potential. Consequently, expanding H\(t) to 0(r) 3 ) around 
x = gives an energy shift e a (t) = fl a (t) sin 2 (0), a linear 
potential of strength g a (t) = rjft a (t) sin(20) and frequency 
change ui a = ujq + Arj 2 Q, a (t) cos(20), where a = {till- 
Accordingly, by choosing the appropriate relative phase, the 
optical dipole potential can cause the oscillator to be tight- 
ened (<f> = 0), slackened (</> = 7r/2) or displaced (<f) = it/ '4), 
while other phases lead to combinations of these effects. For 
concreteness, we focus on a pure displacement quench where 
the perturbation reduces to V — xo(a^ + a) in addition to 
the energy shift e a (t). Since g a (i) oc fl a (t) the protocol can 
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be implemented by varying the laser intensities of the two or- 
thogonally polarized standing waves. 

To examine the verification of the Crooks relation in Eq. (Q3 
we have numerically computed Xf( u ) an d Xs( u ) f° r two dif- 
ferent quenches in a possible 40 Ca + ion experiment 1119 , 



Two experimental limitations were modelled in this calcula- 
tion; First, we used a realistic sampling rate for the mea- 
surement of xf( u ) to account for discrete data. Second, 
an enveloping factor cxp(— u/t), with a decay time r, was 
added to the measurement signal to account for decoherence 
of the entangled state that appears within the scheme I2H1 . 
In Fig. 12 a) the resulting quantum work distributions Pp(W) 
and Pr(W) are shown for a simple tanh forward quench 
from A; = 0, like that shown in Fig. |TJb), described by 
X(t) = A f [l + tanh(i/T)]/2. The switching time-scale T 
used was a fraction of the natural trap frequency 2tt/u>o. Both 
Pf(W) and Pb(W) are composed of 5-peaks, separated by cjo 
and broadened into a continuous spectrum by the decoherence 
envelope. Once inverted, as Pb(— W), the peaks in the two 
spectra line up. As the quench is non-adiabatic, the first-order 
peaks are visible, though much weaker than the carrier peak. 
In Fig. |2|b) the same spectra are plotted for a quench com- 
posed of repeated fast and slow tanh switches (see Fig. Q2c)) 
with the fast switching on the order of a hundredth of 2tt/luo- 
The stronger first-order peaks and now-visible second-order 
peaks evidence how this quench is more non-adiabatic. In 
Fig- I2t c ) we model an additional limitation by computing the 
spectra for the same quench as Fig. 0b) with 0.5% Gaussian 
noise added to Xf( u ) an d Xb( u )- Despite this white-noise the 
first-order peaks remain visible. 

The analysis is completed by extracting the amplitudes 
from Pp(W) and P^{W) of all identifiable peaks. The ra- 
tio Pp(W) / Pb(—W) is then computed for these energies W. 
The Crooks relation in Eq. (Q3 predicts that they will lie on 
a exponential curve, which for the exact (3 and AF of the 
quench is shown on a logarithmic scale in Fig. |2jd) as the 
solid line. The extracted ratios for the three examples given 
in Fig. Ela)-(c) are also plotted and found to cluster tightly 
on this line. For each case, a fitting to exp(AW — B) 
is made with the parameter A providing an estimate of /3, 
thereby establishing that the interferometric protocol also acts 
as a thermometer. This value can be independently com- 
pared to a direct measurement of the initial phonon distri- 
bution of the ion. Using A, the fit parameter B/A subse- 
quently allows an estimate of AF to be extracted lEill . For 
the cases in Fig.|2ja)-(b) the fits essentially yield the exact re- 
sult AF = e a {tR) — ^9<j{tfi) 2 /(jJo, demonstrating the quench 
independence of the relation Eq. ([TJ. The noisy spectra in 
Fig. |2jc) also provides a good estimate (see caption) of both 
/3 and AF from its zeroth- and first-order peaks, illustrating 
a degree of robustness in the scheme when sufficiently non- 
adiabatic quenches are employed. 

Discussions and conclusions - We have outlined a general 
experimental scheme to extract the full statistics of work done 
on a quantum system. Our scheme uses Ramsey interferome- 
try of a single probe qubit to extract the characteristic function 
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FIG. 2: Here a 4c, Ca + is assumed to be confined with an axial trap- 
ping frequency cjq = 300 kHz and in an initial thermal state with 
n = 1. The standing wave optical dipole potential for either po- 
larization are taken as being generated by a 397 nm laser ~ 30 
GHz detuned from the S\i% — P1/2 transition, giving r\ = 0.33, 
and with a maximum Rabi frequency of ft = 150 kHZ fl9|,|20|l. 
A sampling rate of 2 MHz, along with a decoherence time scale 
of r = 15 x (2-k/uq) — 50 /is has been used and measure- 
ments were performed up to a time ~ 500 /is where signal was 
completely damped, (a) The quantum work distributions Pp(W) 
(solid) and Pb(—W) (dashed) are plotted on a logarithmic scale 
for a single tanh switching of A(f), as shown in Fig. [TJb), with 
T — 0.3 x (2n/u}o) = 1 [is. The distributions are continuous due to 
decoherence induced broadening. The dots with horizontal lines de- 
note the peak amplitudes identified from both Pf(W) and Pb ( — W). 
The vertical dashed line is at a frequency of 67 kHz corresponding to 
the exact AF. (b) The plot as (a) but with a repeated tanh switching 
quench (as shown in Fig.[T]c)) with Tf as t = 0.2 x (2tt/ujo) = 0.03 
[is and T s i ow = 3 x (2-k/lvo) = 20 [is. (c) The plot and setup 
as (b) with 0.5% Gaussian noise added to the signals Xf(w) and 
Xb(u). (d) The ratio of the P F (W) / P B (-W) against W evalu- 
ated from the peaks identified in (a) (+), (b) (•) and (c) (x) plot- 
ted on a logarithmic scale. The solid line is Crooks relation from 
Eq. Q} expected from the exact /3 and AF. A best fit of the function 
exp(AW — B) for the noisy spectrum in (c) gives A — 0.72/ojo and 
B/A = 0.20cjo, which compare to the exact values ft = 0.69/ojo 
and AF = 0.22i^o, respectively. 



of the work distribution following an non-equilibrium quench 
of a quantum system. This bypasses the requirement of other 
proposals to implement resource intensive projective energy 
measurements, with the added benefit of being generally ap- 
plicable to a wide range of current quantum technologies. We 
have demonstrated the feasibility of our scheme using a con- 
ventional ion-trap system and standard tools for laser manipu- 
lation under realistic conditions. As such our proposal should 
pave the way for the first experimental verification of fluc- 
tuation relations in the quantum regime. Beyond this, we 
propose that our work is easily adjusted to probe manybody 
systems where recent studies have shown that the statistics of 
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work can be used to shed light on the universal critical fea- 
tures of models from many -body physics Il23l - l25ll . In addition 
further studies have also established an intriguing relationship 
between energy fluctuations in a local quench and the block 
entanglement of a many-body state [2(| and it has been sug- 
gested that these ideas can be tested by a local probe 127 1. 
We therefore suggest that the scheme presented here may be 
used to probe non-local correlations in a macroscopic out-of 
equilibrium many-body system. 
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